Fabrication of Multi-Walled Structure through Parametric Study of Bead Geometries of GMAW-Based WAAM Process of SS309L

In the present study, an attempt is made to investigate and optimize the bead geometries of bead width (BW) and bead height (BH) of SS-309L using an SS316L substrate by employing a gas metal arc welding (GMAW)-based wire-arc additive manufacturing (WAAM) process. The Box–Behnken design approach was used to conduct the trials of single-layer depositions with input variables of travel speed (TS), voltage (V), and gas mixture ratio (GMR). The developed multi-variable regression models were tested for feasibility using ANOVA and residual plots. The data obtained indicated that V had the most significant impact on BW, followed by TS and GMR. For BH, TS had the most significant impact, followed by GMR and V. The results of single-response optimization using a passing vehicle search (PVS) algorithm showed a maximum BH of 9.48 mm and a minimum BW of 5.90 mm. To tackle the contradictory situation, a multi-objective PVS algorithm was employed, which produced non-dominated solutions. A multi-layered structure was successfully fabricated at the optimal parametric settings of TS at 20 mm/s, of voltage at 22 V, and of GMR at 3. For multi-layer structures, fusion among the layers was observed to be good, and they were found to be free from the disbonding of layers. This revealed the suitability of the PVS algorithm for generating suitable optimal WAAM variables. We consider the current work highly beneficial for users fabricating multi-layer structures.


Introduction
Additive manufacturing (AM) is now preferred as a widely accepted technique over traditional manufacturing methods because it produces near-net-shaped components quickly and efficiently using various materials [1,2]. Presenting opportunities while minimizing expenses and enhancing efficiency has become increasingly appealing to manufacturers seeking to optimize their processes [3]. As a result, AM has gained significant traction in recent years, with many companies adopting it on a large scale to gain a competitive edge in their respective industries [4]. The aerospace, automotive, and biomedical areas are becoming more interested in metal additive manufacturing [5]. Wire-arc additive manufacturing (WAAM) is a technique that can save time, costs, and material when making components compared to other AM methods [6,7]. Three key categories used to classify AM for heat sources are electron beams, laser beams, and electric arcs. Electron and laser beams require metal powder as feedstock material, which restricts their production capability [8,9]. Due to this reason, these two techniques have certain limitations for largerscale productions [10]. However, using an electric arc as a heat source is mainly suitable for producing intricate and complex large-scale components at lower cost and reducing material waste due to the more significant deposition rate [11,12]. The electric arc method methodology as compared to the conventional one. Mishra et al. [31] used the optimization approach for simultaneous topology and deposition direction in WAAM using a mathematical model and a combination of GA-and gradient-based optimization techniques. This resulted in an improvement in part quality and a reduction in manufacturing time compared to conventional methods.
Stainless steel (SS) has increased its popularity in WAAM processing. SS is highly valued for its capacity to resist corrosion, rust, and staining [32]. Owing to its excellent strength, durability, and aesthetic appeal, it is primarily preferred for various applications [24]. It is also used in the medical and aerospace industries due to its biocompatibility and excellent ratio of strength-to-weight [33]. SS309L is a heat-resistant alloy commonly used in the chemical and petrochemical industries. SS309L has increased carbon presence compared to other steels, which gives it improved high-temperature strength [34]. The alloy is also resistant to suffixation and carburization, which makes it suitable for use in environments where these processes occur [35][36][37].
Based on past studied work, minimal work has been presented on the experimental investigation and optimization of process parameters for bead morphologies using the GMAW-based WAAM process for SS309L. Thus, the current study attempts an experimental investigation and optimization of the bead geometries of BW and BH of SS-309L using an SS316L substrate by employing a GMAW-based WAAM process. The optimized set of parameters is used to fabricate a multi-walled structure. TS, V, and GMR are elected as WAAM variables based on preliminary experimental trials, machine limits, and recently studied work. At the same time, BH and BW are taken as responses of bead morphology. A BBD is used to generate an experimental matrix for single-layer depositions, and the obtained results are analyzed through ANOVA, residual plots, and main effect plots. WAAM variables are optimized through a PVS algorithm using the empirical relations developed through the BBD. The multi-walled structure is then successfully fabricated at optimal parametric settings. We consider the current work with optimized parametric settings to be highly beneficial for users fabricating multi-layer structures.

Experimental Setup and Plan
In the present study, a metal wire of SS309L with a 1.2 mm diameter was used, and the bead on the plate was deposited on stainless-steel 316L-grade substrate plates using GMAW-based WAAM. Table 1 shows the chemical compositions of the filler wire and substrate, respectively. The GMAW process uses a temporary heat source to heat, melt, and solidify two parent metals and a filler material in a limited fusion zone to make a joint between the parent metals. An autonomous wire feeder constantly feeds the wire electrode via the tip of the torch, where the heat from the welding arc melts it. The distance between the end of the melted electrode and the molten weld pools and the transfer of molten metal to the weld pools controls the heat. The GMAW welding parameters determine the quality and cost of the welded joint. An ideal arc is formed if all the welding parameters are optimal and in accordance. Figure 1 displays the experimental setup used in the present study for the GMAW-based WAAM process.  The setup used in the study had the components of a wire feeder, a GMAW torch, shielding gas, a mixing chamber, and a controller. The build volume for the machine used was 220 × 220 × 500 mm. The torch was enabled to move in the x, y, and z axes to deposit the material on substrate plates. The controller was provided with input through G-code programming via a computer interface in the experimental setup. COLTON iFLEX 350 was used as a power source to heat and melt the metal wire.
A BBD was used to generate an experimental matrix for single-layer depositions on substrate plates using metal wire of SS309L. By using orthogonal arrays of a BBD, multiple factors can be tested with minimal experimental runs [38,39]. This approach helps systematically identify the essential factors that affect product or process quality [40]. A BBD also gives the relationships between a response variable and multiple input variables, such as polynomial regression models, to identify optimal input conditions that maximize the response variable [41,42]. TS, V, and GMR were elected as WAAM variables based on preliminary experimental trials, machine limits, and recently studied work. At the same time, BH and BW were taken as responses of bead morphology. GMR represented the proportions of CO2 gas, and the remainder was argon. Preliminary trials were carried out to identify the range of the selected variables following the BBD. Throughout the single-layer WAAM depositions, a weld bead length of 150 mm, an arc length of 3 mm, and a gas flow rate of 15 L/min were maintained. Table 2 displays the machining conditions used in the present study. The setup used in the study had the components of a wire feeder, a GMAW torch, shielding gas, a mixing chamber, and a controller. The build volume for the machine used was 220 × 220 × 500 mm. The torch was enabled to move in the x, y, and z axes to deposit the material on substrate plates. The controller was provided with input through G-code programming via a computer interface in the experimental setup. COLTON iFLEX 350 was used as a power source to heat and melt the metal wire.
A BBD was used to generate an experimental matrix for single-layer depositions on substrate plates using metal wire of SS309L. By using orthogonal arrays of a BBD, multiple factors can be tested with minimal experimental runs [38,39]. This approach helps systematically identify the essential factors that affect product or process quality [40]. A BBD also gives the relationships between a response variable and multiple input variables, such as polynomial regression models, to identify optimal input conditions that maximize the response variable [41,42]. TS, V, and GMR were elected as WAAM variables based on preliminary experimental trials, machine limits, and recently studied work. At the same time, BH and BW were taken as responses of bead morphology. GMR represented the proportions of CO 2 gas, and the remainder was argon. Preliminary trials were carried out to identify the range of the selected variables following the BBD. Throughout the single-layer WAAM depositions, a weld bead length of 150 mm, an arc length of 3 mm, and a gas flow rate of 15 L/min were maintained. Table 2 displays the machining conditions used in the present study.
The selected input parameters were varied at 3 levels, and 15 runs were carried out following the BBD matrix. The bead geometries of BH and BW were investigated for all the experimental trials. Figure 2 displays the single-layer depositions of fifteen trials by following the BBD matrix, as shown in Table 3 (run order). All the bead-on-plate samples were visually checked and found to be free of any lack of fusion, porosity, or any such defects. This established the workable parameter range of the selected study.  The selected input parameters were varied at 3 levels, and 15 runs were carried out following the BBD matrix. The bead geometries of BH and BW were investigated for all the experimental trials. Figure 2 displays the single-layer depositions of fifteen trials by following the BBD matrix, as shown in Table 3 (run order). All the bead-on-plate samples were visually checked and found to be free of any lack of fusion, porosity, or any such defects. This established the workable parameter range of the selected study.    Optical microscopy was used to measure the bead morphologies of each single-layer deposition. Figure 3 represents the methodology used for the determination of bead geometries. Three different cross-sections of the bead depositions were cut to measure the bead geometries at various locations of the deposition. For better accuracy and more reliable results, their average values were considered in the present study. Optical microscopy was used to measure the bead morphologies of each single-la deposition. Figure 3 represents the methodology used for the determination of bead ometries. Three different cross-sections of the bead depositions were cut to measure bead geometries at various locations of the deposition. For better accuracy and more r able results, their average values were considered in the present study.

Optimization Using PVS Algorithm
Savsani and Savsani [43] studied the passing vehicle search (PVS) algorithm, par ularly for design problems related to engineering. The algorithm imitates the passing o vehicle on a two-lane highway, emphasizing the principle of safe overtaking opportu ties. The mechanism involves several interdependent and complex parameters, includ the availability of gaps in oncoming traffic, the speed and acceleration of individual ve cles, road conditions, overall traffic, weather, and driver skill. While the PVS algorith has demonstrated its usefulness, its application in real-world scenarios with various co plexities and uncertainties still requires further investigation and refinement. The al rithm presented in this paper accounts for three types of vehicles (Oncoming Vehicle OV; Front Vehicle-FV; and Back Vehicle-BV) on a two-lane highway. When a Back V hicle desires to overtake a Front Vehicle, it must have a faster speed than the latter. Ov taking cannot occur if the Back Vehicle has a lower speed than the Front Vehicle. The gorithm also considers the speed and position of the Oncoming Vehicle and its respect distances and velocities when determining the feasibility of an overtaking maneuver.

Optimization Using PVS Algorithm
Savsani and Savsani [43] studied the passing vehicle search (PVS) algorithm, particularly for design problems related to engineering. The algorithm imitates the passing of a vehicle on a two-lane highway, emphasizing the principle of safe overtaking opportunities. The mechanism involves several interdependent and complex parameters, including the availability of gaps in oncoming traffic, the speed and acceleration of individual vehicles, road conditions, overall traffic, weather, and driver skill. While the PVS algorithm has demonstrated its usefulness, its application in real-world scenarios with various complexities and uncertainties still requires further investigation and refinement. The algorithm presented in this paper accounts for three types of vehicles (Oncoming Vehicle-OV; Front Vehicle-FV; and Back Vehicle-BV) on a two-lane highway. When a Back Vehicle desires to overtake a Front Vehicle, it must have a faster speed than the latter. Overtaking cannot occur if the Back Vehicle has a lower speed than the Front Vehicle. The algorithm also considers the speed and position of the Oncoming Vehicle and its respective distances and velocities when determining the feasibility of an overtaking maneuver.
On a two-lane road, there are three different vehicles (OV, FV, and BV) with different velocities (V 1 , V 2 , and V 3 ), and their respective distances can be determined at any given time. The distance between BV and FV can be denoted as 'x', and that between FV and OV can be marked as 'y'. Their velocities impose a primary constraint, where either FV's velocity is slower than BV's (V 1 > V 3 ) or BV's velocity is slower than FV's. If FV exceeds the speed of BV, overtaking becomes impossible, and BV maintains its desired velocity. It is possible to pass only if FV's velocity is less than that of BV. An additional requirement for passing is that the distance between FV and BV during the maneuver is shorter than the distance covered by OV within the same timeframe. Thus, different conditions emerge for the selected vehicles. Employing a human-activity-based technique, the PVS algorithm models the passing behavior of vehicles and offers a meta-heuristic optimization approach capable of finding optimal or near-optimal solutions for given objective functions.

Results and Discussion
The BBD matrix of the RSM approach was utilized to obtain the results of BW and BH, as shown in Table 3. The most suitable condition for fabricating a multi-walled structure is to achieve maximum BH and minimum BW. The study also generated multi-variable non-linear regressions between the machining factors and the responses.

Empirical Relations for BH and BW
Multi-variable non-linear regressions were generated using Minitab v17 with the BBD of RSM to establish the relationships between the WAAM variables and bead geometries (BH and BW). These equations provided a starting point for evaluating response values beyond the experimental matrix of the BBD. Equations (1) and (2) are the regression equations for BH and BW, respectively, which were derived using the stepwise method of statistical approach with Minitab software.

ANOVA for BW and BH
A statistical method, analysis of variance (ANOVA), was used to determine the significance of the factors affecting the response variables. It measured the variability between different levels of the factors to determine which factors significantly impacted the response variables. The adequacy and reliability of the resulting regression equations were tested through an ANOVA analysis. Minitab v17 software was utilized to assess the significance of the model terms at a confidence level of 95%. Terms with probability values less than 0.05 were considered to significantly impact the response variables, while non-significant terms were treated as irrelevant [44]. Table 4 showcases the ANOVA results for BW and BH, and Table 5 shows the results of the model summaries for BW and BH. The statistical analysis of the output factors of BW and BH revealed significant contributions from the regression, linear, square, and two-way interaction models.  The three WAAM variables of V, TS, and GMR were all significant factors in the case of BW, while for the BH response, TS and GMR were substantial variables. TS had the highest contribution to BH, while V was observed to have the highest effect on BW. The small impact of the error term on all the responses indicated high accuracy in predicting values with minimal errors, and lack of fit was statistically non-significant, confirming the accuracy of the ANOVA results [45]. The model was deemed appropriate for predicting the output value as a result. The generated regression equations were deemed reliable and dependable for predicting the values of BW and BH, as evidenced by the significance of the model terms and a close-to-one R 2 value, indicating effective prediction. The model's effectiveness was evaluated by analyzing its R 2 values. The R 2 values obtained for BW and BH were 0.9921 and 0.9578. Such high R 2 values suggest that the model's predictions were accurate and closely matched the actual data. The model was a good fit for the data, as the R 2 values were close to one.

Residual Plots for Output Measures
To ensure the validity of the ANOVA model, specific assumptions must be met, and residual plots must be utilized to confirm the analysis outcomes [45]. Residual plots were used for the validation of the statistical results given by the ANOVA. The residual plots contained four individual analyses. Successful validation of these four plots confirmed the adequacy and reliability of ANOVA results and developed regression equations. Additionally, these validations meant that the regression equations could be used for the future prediction of outcomes for any parametric levels within the range of input parameters. The residual plots for BH, consisting of four plots, are presented in Figure 4.
To ensure the validity of the ANOVA model, specific assumptions must be met, and residual plots must be utilized to confirm the analysis outcomes [45]. Residual plots were used for the validation of the statistical results given by the ANOVA. The residual plots contained four individual analyses. Successful validation of these four plots confirmed the adequacy and reliability of ANOVA results and developed regression equations. Additionally, these validations meant that the regression equations could be used for the future prediction of outcomes for any parametric levels within the range of input parameters. The residual plots for BH, consisting of four plots, are presented in Figure 4. The normality plot exhibited a linear trend that supported the ANOVA model's suitability, indicating that the residuals followed a normal distribution. The versus fit plot also showed that the fits were randomly distributed around the source. In contrast, the histogram plot displayed a bell-shaped curve, which indicated the ANOVA data well. Furthermore, the absence of any particular pattern in the versus order plot further confirmed the ANOVA statistics, leading to better predictions of future outcomes. Figure 5 shows similar findings in the case of BW response. Thus, this shows that the generated regression equations were found to be adequate and reliable for future predictions of outcomes. The normality plot exhibited a linear trend that supported the ANOVA model's suitability, indicating that the residuals followed a normal distribution. The versus fit plot also showed that the fits were randomly distributed around the source. In contrast, the histogram plot displayed a bell-shaped curve, which indicated the ANOVA data well. Furthermore, the absence of any particular pattern in the versus order plot further confirmed the ANOVA statistics, leading to better predictions of future outcomes. Figure 5 shows similar findings in the case of BW response. Thus, this shows that the generated regression equations were found to be adequate and reliable for future predictions of outcomes.

Main Effect Plots for Bead Width and Bead Height
The main effect plot, as shown in Figure 6, demonstrated the trends followed by V, TS, and GMR variance for bead height response, and it can be observed that a slight decrement in the BH of the deposited material was observed with an increase in V.

Main Effect Plots for Bead Width and Bead Height
The main effect plot, as shown in Figure 6, demonstrated the trends followed by V, TS, and GMR variance for bead height response, and it can be observed that a slight decrement in the BH of the deposited material was observed with an increase in V.
This phenomenon could be attributed to the increased heat generated through electrical resistance heating and the absorption of electrons onto the wire tip, providing the necessary energy for the melting and superheating of the wire electrode material [46]. The TS showed a trend wherein its increase led to decreased BH. This was because the higher movement speed of the torch allowed less time for deposition, leading to lesser BH [27,47].
The graph of BH vs. GMR shows that, with an increase in GMR value, the BH increased, and it also became obsolete when talking in relation to BH, as it just safeguarded the pool against air contamination. As per the main effect plots, the WAAM variables of V at 22 V, TS at 16 mm/s, and GMR at 1 were desirable to obtain the maximum BH value.

Main Effect Plots for Bead Width and Bead Height
The main effect plot, as shown in Figure 6, demonstrated the trends followed by V, TS, and GMR variance for bead height response, and it can be observed that a slight decrement in the BH of the deposited material was observed with an increase in V. This phenomenon could be attributed to the increased heat generated through electrical resistance heating and the absorption of electrons onto the wire tip, providing the necessary energy for the melting and superheating of the wire electrode material [46]. The TS showed a trend wherein its increase led to decreased BH. This was because the higher movement speed of the torch allowed less time for deposition, leading to lesser BH [27,47]. The graph of BH vs. GMR shows that, with an increase in GMR value, the BH increased, and it also became obsolete when talking in relation to BH, as it just safeguarded the pool against air contamination. As per the main effect plots, the WAAM variables of V at 22 V, TS at 16 mm/s, and GMR at 1 were desirable to obtain the maximum BH value.
The main effect plot, as shown in Figure 7, demonstrates the trends followed by V, TS, and GMR variance for bead height response. The plot of BW vs. V shows that the BW The main effect plot, as shown in Figure 7, demonstrates the trends followed by V, TS, and GMR variance for bead height response. The plot of BW vs. V shows that the BW of deposition increased with an increase in voltage from 22 V to 26 V. This was due to higher voltage making wide arcs, leading to significant drops of molten metal [48]. The TS trend shows that the BW decreased with an increase in TS. This happened because, with higher TS, less metal could be dropped on the same position, leading to less deposition and, thus, smaller BW [25]. Finally, with higher values of GMR, the BW caused an insignificant change as, at first, it decreased and then increased, negating the changes or uncertainties caused. The gas mixture ratio, on the other hand, had minimal impact on the BW response. Its primary function was to protect the weld pool from atmospheric contamination, and it did not significantly affect the weld BW. As per the main effect plots, the WAAM variables of V at 22 V, TS at 24 mm/s, and GMR at 5 were desirable to obtain the minimum BW value. of deposition increased with an increase in voltage from 22 V to 26 V. This was due to higher voltage making wide arcs, leading to significant drops of molten metal [48]. The TS trend shows that the BW decreased with an increase in TS. This happened because, with higher TS, less metal could be dropped on the same position, leading to less deposition and, thus, smaller BW [25]. Finally, with higher values of GMR, the BW caused an insignificant change as, at first, it decreased and then increased, negating the changes or uncertainties caused. The gas mixture ratio, on the other hand, had minimal impact on the BW response. Its primary function was to protect the weld pool from atmospheric contamination, and it did not significantly affect the weld BW. As per the main effect plots, the WAAM variables of V at 22 V, TS at 24 mm/s, and GMR at 5 were desirable to obtain the minimum BW value.

Optimization of BH and BW Responses using PVS Algorithm
The main effect plots of the response variables demonstrated the contradictory levels of WAAM variables for the optimal levels of BH and BW responses. This meant that the PVS algorithm was employed to optimize the responses. During the execution of the PVS

Optimization of BH and BW Responses Using PVS Algorithm
The main effect plots of the response variables demonstrated the contradictory levels of WAAM variables for the optimal levels of BH and BW responses. This meant that the PVS algorithm was employed to optimize the responses. During the execution of the PVS technique, the upper and lower bounds of the design variables were considered as follows: 22 ≤ V ≤ 26; 1 ≤ GMR ≤ 9; and 16 ≤ TS ≤ 24. Table 6 shows the results of the single-response optimization. Apart from voltage, the other two WAAM variables were observed at different levels. Additionally, the BW response value increased for the maximization of the BH response, which was not desirable, and vice versa. To validate the findings of the PVS algorithm, actual experiments were carried out. The smallest error between the results achieved from the PVS algorithm and validation trials resulted in good agreement between the bead geometry and the WAAM variables. This demonstrated the suitability of the PVS algorithm for the developed regression models.
Although the responses of BW and BH exhibited conflicting tendencies concerning the levels of WAAM design variables, it was crucial to determine the optimal design variable combination that could enhance both responses. This was necessary to achieve the desired quality in both BW and BH simultaneously. Owing to such reasons, multi-objective optimization was necessary. Table 7 shows the results of such an optimization, or Pareto points. Pareto fronts were employed to determine non-dominated solutions that satisfied the requirements of various industrial applications.

Fabrication of Multi-Walled Structure
An objective function was selected to fabricate a multi-layered structure by assigning equal importance to the BH and BW responses. The PVS algorithm was used for obtaining the levels of the WAAM variables, and it showed at the voltage of 22 V, the TS of 20 mm/s, and the DOP of 3 response values of BH at 6.46 mm and BW at 6.33 mm for single-layer deposition. A multi-layered structure was fabricated at these parametric settings. The multi-layer structure was fabricated at the optimal parameter settings of WAAM variables, as shown in Figure 9. For better geometry accuracy, the multi-layer structure was fabricated through layer-on-layer deposition following 180-degree turns of filler wire.

Fabrication of Multi-Walled Structure
An objective function was selected to fabricate a multi-layered structure by assigning equal importance to the BH and BW responses. The PVS algorithm was used for obtaining the levels of the WAAM variables, and it showed at the voltage of 22 V, the TS of 20 mm/s, and the DOP of 3 response values of BH at 6.46 mm and BW at 6.33 mm for single-layer deposition. A multi-layered structure was fabricated at these parametric settings. The multi-layer structure was fabricated at the optimal parameter settings of WAAM variables, as shown in Figure 9. For better geometry accuracy, the multi-layer structure was fabricated through layer-on-layer deposition following 180-degree turns of filler wire. A cooling time of 60 s was applied between successive layers to reduce residual stresses. Fusion among the layers was observed to be good, and it was found to be free from the disbonding of layers. Some extra lumps of metal core were noticed on the extreme sides of the structure. However, these were effectively eliminated in post-processing. This revealed the suitability of the PVS algorithm for generating the suitable optimal WAAM variables. Therefore, the present work effectively demonstrated the requirement of having optimal parametric settings and the necessity of parametric optimization for fabricating thin, multi-walled structures using a GMWA-based WAAM process for SS-309L using an SS316L substrate. We believe that the present work may be useful for researchers and industrial applications to find optimal sets of parameters.

Conclusions
The present study used a GMAW-based-WAAM process for SS309L wire. V, TS, and GMR were identified as machining parameters, and BH and BW were output factors. Experiments were conducted following a Box-Behnken design. Optimization of bead geometries was obtained through the application of a PVS algorithm. Based on the key findings and results, the following conclusions can be drawn: • Multi-variable non-linear regressions were generated among the WAAM variables and output responses. • ANOVA was employed to validate the appropriateness and reliability of the obtained regression equations. The ANOVA revealed that the quadratic model, including linear, squared, and interaction model terms, was statistically significant for both the bead height and width responses. The lack of fit results signified the model's suitability and acceptability for both responses. The validation results from the ANOVA of an R-squared value close to one showed that the model was adequate and acceptable.
• TS had the largest impact on BH response, followed by GMR, while TS followed V and GMR showed a substantial impact on BW response.

•
The influences of the WAAM variables (TS, V, and GMR) were studied on the BH and BW responses. It showed conflicting situations for attaining the desired levels of bead geometries. A cooling time of 60 s was applied between successive layers to reduce residual stresses. Fusion among the layers was observed to be good, and it was found to be free from the disbonding of layers. Some extra lumps of metal core were noticed on the extreme sides of the structure. However, these were effectively eliminated in post-processing. This revealed the suitability of the PVS algorithm for generating the suitable optimal WAAM variables. Therefore, the present work effectively demonstrated the requirement of having optimal parametric settings and the necessity of parametric optimization for fabricating thin, multi-walled structures using a GMWA-based WAAM process for SS-309L using an SS316L substrate. We believe that the present work may be useful for researchers and industrial applications to find optimal sets of parameters.

Conclusions
The present study used a GMAW-based-WAAM process for SS309L wire. V, TS, and GMR were identified as machining parameters, and BH and BW were output factors. Experiments were conducted following a Box-Behnken design. Optimization of bead geometries was obtained through the application of a PVS algorithm. Based on the key findings and results, the following conclusions can be drawn:

•
Multi-variable non-linear regressions were generated among the WAAM variables and output responses. • ANOVA was employed to validate the appropriateness and reliability of the obtained regression equations. The ANOVA revealed that the quadratic model, including linear, squared, and interaction model terms, was statistically significant for both the bead height and width responses. The lack of fit results signified the model's suitability and acceptability for both responses. The validation results from the ANOVA of an R-squared value close to one showed that the model was adequate and acceptable. • TS had the largest impact on BH response, followed by GMR, while TS followed V and GMR showed a substantial impact on BW response.

•
The influences of the WAAM variables (TS, V, and GMR) were studied on the BH and BW responses. It showed conflicting situations for attaining the desired levels of bead geometries. • Single-response optimization using the PVS technique obtained a maximum BH and a minimum BW of 9.48 mm and 5.90 mm, respectively. Pareto fronts were employed to determine non-dominated solutions that satisfied the requirements of various industrial applications.

•
The multi-layered structure was successfully fabricated at optimal parametric settings of V at 22 V, TS at 20 mm/s, and DOP at 3. For the multi-layer structure, fusion among the layers was observed to be good, and it was found to be free from the disbonding of layers. This revealed the suitability of PVS for generating suitable, optimal WAAM variables.

•
The present work effectively demonstrated the requirement of having optimal parametric settings and the necessity of parametric optimization for the fabrication of thin, multi-walled structures using a GMWA-based WAAM process for SS-309L using an SS316L substrate. We believe that the present work may be useful to researchers' industrial applications to find optimal sets of parameters. As a result, in future work, fabricated multi-layered structures will be examined using microstructure investigations and mechanical properties such as tensile, impact, and microhardness testing.